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Abstract. We continue our study of the quantum mechanical motion in the x^y^ 
potentials for n = 2,3, which arise in the spatially homogeneous limit of the Yang-Mills 
(YM) equations. In the present paper, we develop a new approach to the calculation of 
the partition function Z{t) beyond the Thomas- Fermi (TP) approximation by adding 
a harmonic (Higgs) potential and taking the limit w — > 0, where v is the vacuum 
expectation value of the Higgs field. Using the Wigner-Kirkwood method to calculate 
higher-order corrections in ?i, we show that the limit w — > leads to power-like 
singularities of the type which reflect the possibility of escape of the particle 

along the channels in the classical limit. We show how these singularities can be 
eliminated by taking into account the quantum fluctuations dictated by the form of 
the potential. 
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1. Introduction 

We here continue our study of tlie quantum mechanical motion in the x'^y'^ potentials 
of phase space dimensions 2n with n = 2,3, which arise in the spatially homogeneous 
limit of the Yang-Mills (YM) equations. As is well known [ij (see p] for a review), these 
systems exhibit a rich chaotic behavior despite their extreme simplicity. Especially the 
n = 2 model, the central object of this and also our previous investigation (see [Hj - we 
will henceforth refer to this work as I), has been widely studied. 

Classically, this model possesses a logarithmically divergent volume of energetically 
accessible phase space Te 0111 t, but its quantum mechanical version (YM quantum 
mechanics, YMQM) has a discrete spectrum ^ [7]. Physically, the explanation is 
obvious: Quantum fluctuations, e.g. zero-point fluctuations, forbid that the trajector 
escapes along the x or y axis where the potential energy vanishes. The system is 
thus confined to a finite volume, and this implies the discreteness of the energy levels. 
Classically this escape is always possible without increasing energy. As we shall see 
below, these classically allowed configurations result in singularities of the quasiclassical 
partition function. 

In I we calculated the higher-order quantum corrections to the partition function 
(heat kernel) Z{t) for the x'^y'^ potential using the approximation jHJ based on 
the adiabatic separation of the motion in x and y in the hyperbolic channels of the 
equipotential surface xy = const. The main assumption of this method is that the 
final results of the calculations do not depend on the artificial boundary Q dividing 
the central region x,y E [—Q,Q] from the channels x,y E [Q,oo]. We showed in I 
that this assumption, after improvement of the quantum mechanical treatment of the 
motion in the channels, is correct not only for the Thomas-Fermi (TF) term but also for 
the leading (in powers of tQ"^ ^ 1) higher-order quantum corrections, and we derived 
Q-independent asymptotic series in the parameter = g'^hH^ for contribution of each 
region to Z{t). 

In the present paper, like in I, we explore the properties of the x'^y'^ model beyond 
the TF approximation, but we pursue a different approach. We calculate the higher- 
order corrections to Z{t) as in ^U] by starting from the Yang-Mills-Higgs quantum 
mechanics (YMHQM) and taking the limit v ^ 0, where v is the vacuum expectation 
value of the Higgs field. In the n = 2, 3 x^y^ models, v determines the strength of the 
harmonic potential 



Due to the above mentioned logarithmic infinity of the classical phase space volume 
of the x^y^ model it is impossible to completely disentangle the nonlinear coupled 

I This is in violation of Weil's famous theorem [1] , which states that the average energy level density 
dN/dE is asymptotically proportional to F^;. 
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oscillators from the harmonic oscillations generated by the Higgs potential leading to 
a term in Z{t) proportional to Inf. In the n = 3 case, which has a finite phase space 
volume at fixed energy, this method yields an expression for Z{t) that coincides in the 
limit f — ^ with the one obtained adiabatic separation method p. This is due to the 
negligible time spent by the classical trajectory in the depth of the hyperbolic channels 
[H]. Higher-order corrections change this situation essentially, as we shall see below. 

Here we use the approach of ref. ^U] with the limiting procedure v ^ for the 
calculation of Z{t) beyond the TF approximation by applying the Wigner-Kirkwood 
(WK) method pTj IT^ IT^ (see |T3] for a review of the WK method). The higher-order 
corrections 0{h^) in the WK approach lead to a new phenomenon in the limit f — > 0: 
for k > 2 they yield power-like singularities of the form f These singularities are not 
cancelled at a given power k as one might expect. The situation is completely different 
when one includes the quantum fluctuations in the channels of the x^y^ potential to 
all orders. These generate a conflning potential, which does not disappear in the limit 
f — >■ 0, closes the flat direction, and eliminates the mentioned singularities, which are 
essentially classical. Taking these fluctuations into account, we are able to compare the 
expression for Z{t) obtained by our method with the result for Z{t) obtained by the 
method of |Hl El in the TF approximation. 

Concerning the n = 3 case with its flnite phase-space volume, the singularities 
appear also at the higher-order corrections in contrast to the TF approximation and 
again are eliminated by the quantum fluctuations corresponding to the speciflc quartic 
form of the potential characteristic of the YM quantum mechanics. § 

Finally, we develop a novel approximation scheme, which relies on a resummation of 
certain terms in the KW expansion and thus introduces a nonvanishing Higgs potential, 
which avoids the divergences of the TF approximation. This approach is motivated by 
the need to take the quantum fluctuations inside the hyperbolic channels into account 
even in the lowest order approximation. We show that the new approach reproduces the 
TF result obtained with the method of jHl El without requiring an artiflcial subdivision 
of the phase space into different regions. A modifled WK expansion can be derived to 
systematically improve on this result. 

In the next two sections we present the YMH system and the WK method of 
calculating Z{t) beyond the TF approximation. 

2. Yang-Mills-Higgs classical and quantum mechanics: The Thomas-Fermi 
approximation 

There are several mechanisms that can suppress and even eliminate the classical chaos 
of the YM system (see |2]). One is the Higgs mechanism jTHj. For spatially homogeneous 
flelds (long wave length limit of YM system), if only the interaction of the gauge flelds 

§ In the corresponding supersymmetric Yang-Mills quantum mechanical system, the wave function 
is not confined due to the cancellations between bosonic and fermionic degrees of freedom, and the 
spectrum is continuous ]15| . 
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with the Higgs vacuum is considered, the classical Hamiltonian for n = 2 is given by 
the expression 



H = l{pl + pl) + 'jxV + j{x' + y% (1) 



where v = (0) is the vacuum expectation value of the Higgs field (p. It is known that 
there is a classical transition from chaos to regular motion as v gets large enough. More 
precisely, the chaos disappears when g'^v^/E > 0.6 ^Hl; where E = H is the energy. 
The analogous transition in the adjacent energy- level spacing distribution was predicted 
[T7] and established in several papers [IB]. The quantized counterpart of (JH) is 



^^ = -TrVL + v^V + v(^' + r). (2) 



As in I, we measure all quantities in units of the energy E with dimensions [H] = 1, [t] = 
-1, [x], [y] = 1/4, [g] = 0, [v] = 1/4, [U] = 3/4. 

It is obvious that the operator has a discrete spectrum as it has for w = 0. 
The TF approximation to the heat kernel or partition function Z{t) = Tr[exp(— ti^)] 
is the standard lowest-order semiclassical approximation valid for small ht^^^ <^ 1. 
It is obtained by substituting the classical Hamiltonian for its quantum counterpart 
and replacing the trace of the heat kernel by the integral over the phase-space volume 
normalized by (27?^)"", where 2n is the phase-space dimension. In other words, the 
TF approximation takes into account only the discreteness of the quantum mechanical 
phase space, but considers momenta and coordinates (in our case, the field amplitudes 
X and y) as commuting variables. This method was used in numerous papers (see 
e.g. 13 EH!). For the calculation of the energy level density p{E) = dN{E) /dE at 
asymptotic energies, the TF approximation is a consistent approach since, as we shall see 
below, all corrections to the TF term are structures with factors h^t^ with k, i positive 
integers. For the asymptotic energy level density p{E) or N{E) these corrections are 
negligible according the Karamata-Tauberian theorems jBJ |7] relating the most singular 
part of Z(t) to the asymptotic level density, N{E) = J dEp{E) = C^^[Z{t)/t], where 
C^^ denotes the inverse Laplace transform. 

In PU] Z{t) and N{E) were calculated for the Hamiltonian (^J). We give the precise 
expression for Z{t) of the YMHQM system in the TF approximation: 

m= ^,L,. expfgUof^V (3) 



/2TTghH^/^ 

where Kq{z) is the modified Bessel function of the third kind. For the most interesting 
limit V ^ we get: 

Zit) - fin ^i^-c], (4) 



^ghh^/^ 

where C is the Euler constant. The impossibility of disentangling the coupled oscillators 
from the uncoupled ones is expressed by the logarithmic divergence of the phase space 
volume for n = 2 as we already mentioned in the Introduction. Below we compare (jH) 
with the corresponding expression for Z{t) obtained in [HI El for the pure x'^y'^ model. 



Coupled Yang-Mills Oscillators II 



5 



Because we shall often encounter the pre-factor appearing in Q and (@]) in the following, 
we introduce the special symbol K for it: 

K = {27TgWy'/^ = {2nXY'^\ (5) 



3. Beyond the TF approximation: The Kirkwood-Wigner expansion 

In the present paper, unhke in I, we apply the WK expansion in all of phase space, 
avoiding the division of the phase space into a central region and hyperbolic channels. 
As we will see below, this poses no problems as long as f 7^ 0. However, singularities 
appear in the limit v 0, unless the WK expansion is modified to include quantum 
fluctuations in a nonperturbative way. 

Since we described the WK method in detail in I, we only give a very brief outline 
here. We start from the equation (1-12), a set of recurrent differential equations for the 
kernels Wk of the partition function Zk{t) at the fc-th order in h: 



Zk{t) 



{2-nny 



dTWu{rJ;t) e 



-tH 



(6) 



where dV = dxdydpxdpy and the classical Hamiltonian H given by (^. We begin with 
the partition function at second order in ^ (A; = 2). Integrating over p^ and Py and 
making use of the symmetry of the Hamiltonian (P) with respect to the interchange 
X ^ y, we obtain: 



dV W, e-*^ = - 



(7) 



222 

- X y 



(8) 



where we have used the notation (for m > n): \\ 

= 4 / dx dy exp -- (v\x^ + y^) + g 

Jo Jo L 2 ^ 

Note that the factors from the KW expansion and from the normalization of the 
phase space volume element have canceled, making Z2 independent of h. Integrating 
over y and introducing the new variable u = g'^x'^/v'^, we obtain: 



27r(2n - roc 



Jo 



with z = tv^/2g'^. The integral over u is related to the Whittaker function Wx^^{z) (see 
HH], equation 9.222.1): 

v^(2n - l)!!(i;2)™-" /2m + l\ . ^ 

Iran = ^ , \ ' e^'^Z {z) 10 

^n+2(^2m+l V 2 / 2^2^'^' 

For m = n the Whittaker function has only a logarithmic divergence at ^ = 0. For 
m — n > 1 power-like singularities appear. For completeness we give Z2{t) explicitly 

Z^it) = !^ , [2 f-l + ^1 .-F f ^ W . . (.) 



12V27r^ti 

Note that the case m = n needs to be calculated separately from the case m > n\ see below. 
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As is easily seen from (fTT|) . ^2(^) has a singularity of the form f ~^ at f = 0. Using the 
limit of the Whittaker function W\^fj_{z) for small we find: 

Ut) - 'K^^^ {V - 0) (12) 

For completeness we give here the Z2 from the KW method in the limit g = 0, i. e. for 
two free harmonic oscillators using the asymptotic form of the Whittaker function: 

Z2{t) = (13) 
Together with the expression for the TF term ITU] we have 

z,{t) + z,it) = ^ (1 - ^nw) , (14) 

which are the first two terms in the Taylor expansion of the exact partition function for 
the two-dimensional harmonic oscillator at small hvt: 



2 sinh -Tivt 
2 



(15) 



4. Higher-order corrections and the limit f ^ for YMHQM 



At higher order h!^ {k > 2) the mathematical structure of terms Wk and changes 
essentially. In the expression (jH} higher powers m and n appear and the difference 
between them increases {m — n > 1), causing the second index of the Whittaker function 
to exceed fi = 1/2. As a result, in the limit z = tv^ /2g'^ the finite sum of the power- 
like singular terms of W\^^{z) begins to play a crucial role. 

A systematic analysis of the higher-order corrections using Mathematica leads to 
the conclusion that there is a correlation between the power of h (for k > 2) and the 
most singular terms in pO|) : m — n = k/2 (due to the symmetry against exchange x ^ y 
it is always possible to put m > n). The case m = n requires special consideration and, 
as found, leads only to logarithmic singularities. 

It is easy to see that the next, less singular terms correspond to the case m — n = 
{k/2) — 21 {£ = 1,2, . . < k/A). For the general expression of Z^'^\t) we need to 
determine the powers of g'^ and t. For g'^ it is simply (7^™. To determine the power of t 
we note that for each k there are always terms without factors of px and Py and minimal 
power of t at given k. For such terms the factor is t^m-n-i the terms containing 
factors of px and Py the power of t does not change after integration over the momenta. 
For the less singular terms with m — n = ^k — 2£ the corresponding factor is i2m-n-i-3^_ 

Now we are in position to write the general expression for the most singular 
terms contributing to the partition function at the order of h^{k = 2,4,6, . . .). After 
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integration over p^^py^x and keeping only terms having m — n = k/2 {k/2 < m < 
k,0 < n < k/2) we obtain: 

Zt'''\t) = K f^^V^'r f^") {2m- k- 1)!!. (16) 



^ vH ^ 

We see that there are strong singularities at f = 0, and further analysis shows that these 
are not cancelled by the summation of all most singular terms at a given k. For the less 
singular terms with m — n = k/2 — 2i with i < k/A we have: 

Due to the appearance of (m — n) in the argument of the gamma function in (|T7j) . the 
case with m = n must be considered separately. It is obvious that these terms have no 
power-like divergences, but only logarithmic singularities like the TF term. We obtain 



ln^-2C-^fm+- 



(18) 



where il){x) is the logarithmic derivative of the Gamma function. Let us briefly discuss 
these results. The power-like singularities in the KW approach for k > 2 are related 
to the possibility of escaping classically along the axis x = or ?/ = where, in the 
limit f — *• 0, the potential energy vanishes. Non-zero v forbids such escape to infinity. 
These singularities affect any classically calculated distribution function, in particular, 
the heat kernel Z{t). They also show that the trajectories lie deep inside the channels 
most of the time. Quantum mechanical fluctuations forbid any escapes along the axes. 

The singularities have more resemblance with the infrared singularities, they 
are related to the behavior at long distances and different from the usual ultraviolet 
divergences connected with the asymptotic expansion in powers of h^. The absence 
of power-like singularities for m = n is easily explained: for m = n the configurations 
dominate along the diagonals (|x| = \y\), whereas for m ^ n the configurations populate 
the channels and have a trend to escape if quantum fluctuations do not forbid this. It 
is thus clear that we have to take into account the effect of the quantum fluctuations on 
the motion inside the channels, which the perturbative expansion of the WK method 
fails to do. 



5. Quantum fluctuations and power-like singularities 

In this Section we will attempt to include quantum fluctuations created by the form of 
x^y^ potential in a framework, which does not rely directly on the adiabatic separation 
of the motion inside the hyperbolic channels as it was elaborated in detail in I. Let 
us consider the motion along the x-axis. The heat kernel for the x'^y'^ potential 
generates a mean spread in y at the position x of the order of 5y ~ {g'^x'^t/2)~^^'^). 
Quantum mechanics dictates that the spread of the conjugate momentum py is at least 
Spy > hg\x\{t/2Y^'^. Analogous relations hold between y and p^- 5px > %|?/| (t/2)^/^ 
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for the motion along the y-axis. We now propose a modification of the WK formahsm, 
which incorporates these relationships into the generating functional for the expansion 
in powers of h. 

To achieve this, we resum the term —^h^t(AV) in the WK operator (I-IO) to all 
orders by writing 

W{r,p;t) = W{r,p;t)exp (-^t'^Av\ . 



(19) 



Inserting this definition into the differential equation (I-IO) for the WK kernel W yields 
an equation for the phase space function W: 



dW 



;2 r 



A + t\VV'f - ^{p- VV) - ^t\AAV) 



+ -(ip- ntvv) ■ V 

Ti 



t^V{AV) ■ V 



—t^iip- MVV) ■ V(A\/) 



W. 



The term in the exponent of (fT^ : 



-^fAV = -^gH\x' + y'] 



(20) 



(21) 



acts like a Higgs potential with v'^g = What is unusual about this term is 

that the effective potential itself is time-dependent. The connection to the argument 
given at the beginning of this section becomes evident, when one recognizes that the 
exponent represents the lower bound associated with the kinetic energy demanded by 
the uncertainty relation: 



l{Spl + 5pl)>^tAV. 



(22) 



It is straightforward to derive a recursion relation analogous to (1-12) for the 
coefficients of the expansion of W in powers of h: % 

dWu 1 



dt 



ip-[w - tiyv)] Wk-i + - A + t'iyvy - 2tv\/ ■ v Wu- 



'\p-W{AV) 



+ 



2 



-VV ■ V{AV) - -V(AV) • V + -AAV 



k-i- 



(23) 



Before we apply this trick to the quartic YM oscillator, it is useful to briefly explore 
how it works for the harmonic oscillator. It is easy to see that ()23|) yields the correct 



expansion for the partition function for the potential V = ^x^. Indeed, by calculating 

% Note that the expansion of W in powers of h does not yield an expansion of Z{t) strictly in powers 
of h because of the nonpolynomial factor containing h in H19|l . 
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Wq, W2, W4 and integrating over p and x, we obtain 

1 

hvt 

5 



24 

5760^ ^ ^ 



(24) 



giving the expansion in powers of {hvt) (up to of the exact partition function of the 
harmonic oscillator, Z{t) = (2 sinh^ft/2)~^. Note that, in this case, Zj^ contributes to 
all powers in the h expansion of Z{t) beginning with h''~^. 

Now consider the non-linear potential ^g^x^y^. The lowest term is easily obtained 
from by substituting v Vcs'- 



Kexp{—gV\Ko{—gH 



2+3 



16 



16- 



K 



In 



1 



9 



+ 51n2-C 



(25) 



which coincides (with logarithmic precision) with the result obtained in [Hlin] for Z{t) 
The second-order term is 

^ g^t 



247r 



{Aho + gHhi] 



with 



dx 



dyx^V^exp 



1 h'^ ^ 

^222. "' 2/ 2 I 2\.2 

-1^9 X y t - —g [x + y )t 



(26) 



(27) 



These integrals correspond to those defined in (jHI) with the substitution t>^ — ^ h^gH/2. 
The exact analytical expression for these integrals in terms of Whittaker functions was 
given in ()10|) . In the following we retain only the first term in the (finite) power series 
expansion of the Whittaker function; later we shall correct for this simplification. Using 
the condition z = g'^h'^t^/8 = <^ 1, allowing us to neglect terms involving InA^ 
compared with terms of the form A~^, we obtain 

~ a/2^ 4 



(10 



'21 



(gHy/^ (hgtr 
4 



(gHf/^ {hgtr 



yielding the result 



Z2{t) = -K. 



(28) 



(29) 



For Z4 we retain only terms with m — n = 2 as the most singular ones in the limit 
f — s> 0, based on the usual arguments. This gives: 

{hgtr 



z^it) 



1 ^ g't- 

— ion ~r J■'\^ 

30 180 



(30) 
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For the needed integrals we obtain in the same approximation (A^ ^1): 

16 



10 



I20 

hi 



16 
48 



{gHfl^ihgty 
Substituting (jHTjl into (jHUj) we obtain: 

127 

Z^{t) = —-K. 



180 

Collecting all results, we have: 



^0+2+4 (i) = K 



, 1 , ^ 427' 

ln-rT^ + 51n2-C + 



(31) 



(32) 



(33) 



g^nH^ ' " ' 180_ 

We improve upon the approximation made above to the Whittaker function and 
consider the contribution of all singular terms in the asymptotic expansion of W^^fiiz) 
for small z: 



m— n — 1 / -j^ 

T(m — n — p)T in + - + p 
p=0 ^ ^ 



A 



2\ P 



(34) 



where again A^ = g'^hH^ . The complete expression for the contribution to Zk {k 
2, 4, 6, . . .) from the most singular terms is 



K 



2'=(2n- 1)!! t^^„ /I 



p=0 ^ 



T{n+\) 

For the less singular terms {m — n = — 2i, i < ^k) we get 

^2fc-4£(2n- 1)!! 



n' + -+ p 



A 



2\ P 



(35) 



r(n + i) 

g r(-i-2/-p)r(„ + -+p) 



,2\ P 



and for the logarithmic term (m = n) we obtain 



7{m=n) 



KX2 



-lnA2 + 31n2-2C-?/; m + 



1 



(36) 



(37) 



It is clear that these contributions again generate an asymptotic series in the small 
parameter A^ = g'^fi^t^ of the following form: 

'=/2 Ok(On — Till 

Z{t) = K -lnA2 + 51n2-C+ 2^ E 4'^'^ 

fc=2,4... n=0 



p=0 



E r -A;-p rU + - + p 



1 



A 



2\ P 



(3^ 
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where the a,[f/^) (n = 0, 1, . . . , |) are the coefficients of the structures {g'^t)^Imn 
(m — n = |) in the expression / dT M4e~^* {k = 2, 4, 6, . . .), analogous to the coefficients 
a^-* (m = 0, 1, . . . , 3n) in (1-71). For /c = 4, e.g., these numbers are: 

^° -30' -180' - 576- ^^^^ 

We note that this asymptotic series closely resembles the one derived in I using a 
completely different approach, splitting the x — y plane into two integration region 
and treating the quantum fluctuations exactly in the region containing the hyperbolic 
channels. The leading logarithmic term is identical in both cases, but it is not clear 
that the constant coincides. The expansion parameter is the same for both series. 
Unfortunately, our inability to find a simple general algorithm for the coefficients, a^^^^^ 
here and in I, has prevented us to compare the two results in detail. The less singular 
terms also lead, after summation over k, to an asymptotic series in the parameter A^. 
The same is true for the terms zj^'"^^ with m = n. 

At the end of Section 4 we commented on the analogy between the limit f — > and 
the infrared problem. We now can make this analogy more precise. The infrared limit 
corresponds to the behavior of the system at large distances or deep inside the channels, 
in the terminology of the paper I. When we compare the limiting forms for f — of 
the quantities /^n derived here with the expression (1-25) for the analogous quantities 
derived for large distances {tQ'^ ^ 1) in I, we obtain the correspondence: 
n2(m-n) _ _rn-_n 1 

Evidently the limit v ^ corresponds to the limit Q —* oo, supporting our claim that 
it represents the infrared limit. 



6. The three dimensional YMHQM 

Finally, we briefly consider the n = 3 case of the YMHQM model. The Hamiltonian for 
n = 3 YMH classical mechanics is: 

H=^ipl+pl+pl) + ^(xV + + z'x'). (41) 

In its quantum counterpart, is replaced with — /i^V^. As we know the contribution 
to Z{t) from the channels is negligible in the TF approximation j^j if we apply the 
condition tQ^ ^ 1. This is due to the fact that deep in a channel, e. g. along the 
X-axis, the energetically accessible phase space volume gets pinched as and not as 
as it for the n = 2 Hamiltonian. This implies that the limit f — >^ is smooth and 
the expressions for Z{t) from and [T^ coincide at f — 0. Already for the second 
correction to the TF term this is not true as we shall see below. In the WK approach 
for Z2 using (1-12) and © we have 

Z2{t) = , ,3 I dpxdpydpzdxdydzW2{p,x;t) e~^\ (42) 
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with W2 from (1-19). Integrating over the momenta, using the symmetry of the potential 
energy in ()41|) . we have 



^1/2 



-9 h + —h 



where integrals Ji and I2 over x, z are given by 



h 



dxdydz x^e 



-V {x,y,z)t , 



roo 

I2 = dxdydz x^iy"^ + z^^e' 
Jo 



with 



V{x,y,z) 



9 



.2^2 



■V{x,y,z)t. 



^2 2\ 



22, 22, 2 2\|^/2| 2, 2^ 

j{X y +y z + +y +z) 



(43) 

(44) 
(45) 

(46) 



Integrating over x and using polar coordinates r, 4> for the integrations over y and z we 
obtain: 



1 /27r\3/2 



A = - 



27r\3/2 



r(ir 



exp (^—j^tg'^r'^ — ^tv'^r'^ 



2^4 



(^2 + (,2^2)3/2 



oo exp 
r dr 



¥9'r' 



16 



(47) 



(4^ 



{v^ -\- ^2j-2^3/2 

where /o(^) denotes the modified Bessel function of the first kind. We see that I\ is 
divergent at r = if f = in contrast to the TF term. 

Again, we obtain a smooth transition in the limit f — > if we make the substitution 
used in Section V before taking the limit f ^ 0. As for the n = 2 model, here also the 
zero-point quantum fluctuations in the channels generate an effective Higgs potential 
\h^g'^t{x^ + y"^ + z"^) and render the integral /i convergent. After this substitution, 
introducing the new variable 

Z (t) - \ 
with A = gt^/'^h^{<^ 1) and 



= Tei9'^r^, we get: 
Jo(A)+4J2(A)] 



(49) 



u du 



(50) 



(A + 8m)3/2- 

Expanding the exponential function in the small parameter A, the integral (I^Uj) can be 
expressed as a finite sum of the generalized hypergeometric functions. Retaining the 
main terms we have 

Jo(A) ^ 



J2W 



4v^ 

r 



32v^' 



Thus, the second-order correction Z2{t) for the n = ?> YM model becomes 



Z2{t) ^ L 



25/4^3/2 



{gwy/' 



(51) 



(52) 
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where 




is the TF term found in CHI • We see from ()52|) that the quantum corrections at the 
order are parametrically enhanced due to the quantum fluctuations generated by the 
effective Higgs potential ^h'^g'^{x'^ + y"^ + z^)t- 

7. Conclusions 

We have shown in I and here that the richness of the classical YM mechanics with 
the x^y'^ potential translates into, and even gets amplified by, the quantum mechanical 
properties of the system. The YM quantum mechanics exhibits a confinement property, 
which strongly influences the quantum mechanical motion in the x^y"^ potential. At 
higher order in h (up to fi^) this results in the vanishing of the leading quantum 
corrections, when we correctly take into account this property for the motion in the 
hyperbolic channels. We convinced that this property survives at higher orders, although 
we did not explicitly demonstrate it. 

Here we calculated the quantum corrections to the partition function by adding 
a Higgs term to the potential, and found that power-like singularities arise in the 
limit f — ^> 0. We associate these essentially classical singularities with the fact that 
the Wigner-Kirkwood expansion does not take into account the effect of the quantum 
fluctuations on the motion within the channels. When these fluctuations, which are 
dictated by the uncertainty relation and the hyperbolic form of the channels, are taken 
into account, the escape along the coordinate axes, which is classically allowed, is 
prohibited, and the singularities disappear. As a result, the Thomas- Fermi term for 
the partition function acquires a renormalization expressed in terms of an asymptotic 
series in the parameter = g'^h'^t^ within both approaches. 

We hope that the lessons we elicited from the present study of the higher-order 
quantum corrections to the homogeneous limit of the Yang-Mills equations will be useful 
for an improved understanding of the internal dynamics of the Yang-Mills quantum field 
theory. 
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